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The single-particle spectral function and the density response of a two band Emery model for 
CuO chains is calculated for large on-site Cu repulsion U and large on-site energy difference A. For 
U ^ U — A '^ t the eigenfunctions are products of charge and spin parts, which allows analytical 
calculation of spectral functions in that limit. For other parameters numerical diagonalization is 
used. The low energy hole carriers are shown to be the one-dimensional analogs of the Zhang- 
Rice singlets. The validity of the one-band model is discussed. The results are relevant to the 
interpretation of photoemission and EELS experiments on SrCu02 and Sr2Cu03 . 
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I. INTRODUCTION 

One of the most intriguing phenomena in one- 
dimensional electron systems is the so called spin-charge 
separation: the low energy excitations are decoupled col- 
lective modes of charge and spin character, which may 
have different velocities, and are referred to as holons 
and spinous, respectively. As a consequence, the spin 
and charge of an added electron will be spatially sepa- 
rated after some time and there are no Fermi-liquid like 
quasiparticles. While the decoupling exists also in the 
weak coupling limitEJ, it is perhaps best understood for 
the strong coupling limit of the Hubbard model, where 
the Bethe ansatz solution tells us that the wave func- 
tions are factorized into a part describing free spinless 
fcrmionsjjcpresenting the charges and a part representing 
the spina3. This allowed the calculation of tlte dynamical- 
spectral functions of the Hubbard model atn and awayQ 
from half filling with excellent resolution. These calcula- 
tions provided an explanation of the origin of the different 
features in the spectral function. 

The most direct test of the theory is to look at the pho- 
toemission spectra of highly anisotropic materials. The 
nearly ideally one-dimensional CuO chainsu in the charge 
transfer insulators SrCu02 and Sr2Cu03 are perfect can- 
didates, given that the typical energy scale for spin and 
charge excitations is large compared to the experimen- 
tal resolution, making the observation of the low energy 
spectrum possible. The absence of bands would indicate 
that we are not dealing with the usual quasiparticles of 
Fermi liquid theory. On the other hand, there are very 
clear theoretical predictions for the photoemission spec- 
trum of a system where spin-charge separation exists, arid 
indeed, recent phottieHiission experiments on SijC|U02lI u 
and photoemissionLHl^ and electron-energy-lossliil exper- 
iments on Sr2Cu03 seem to indicate that the dynamics 
at low energies can be understood within an effective t-J 
or Hubbard model. 

In comparing the measured spectra with the theoret- 
ical ones, we face the following difficulties: i) the ac- 



tual material is a charge transfer insulator, while the 
Hubbard/i-J model is a Mott insulator. Therefore one 
is lead to question how much of the spectra can be at- 
tributed to generic features where the details of the model 
are not important; ii) For the Cu02 plane, the t-J model 
is derived to describe t 
the Zhang- Rice singlets 



dynamics of complex objects - 
^^' In the CuOs chains the O ions 
are not all identical, and the original picture of Zhang 
and Rice has to be refined; iii) On the theoretical side, 
apart from numerical calculations of the dynamical cor- 
relations which are difficult to interpret, exact and/or 
analyticaL-ccsults are very rare concerning the spectral 
functionBEj^Ej and the optical conductivity (e.g. for the 
Hubbard model see Ref. nfl). 
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FIG. 1. The CuOa chain with 3d orbitals of Cu and 2p 
orbitals of O. The side oxygens (dashed lines) are omitted in 
the two-band model. 

To answer these questions, we will consider the ID 
model involving the Cu ions and the oxygens be- 
tweeur-thetn (see Fig. P, the so called two-band Emery 
mode]ll3't2l, as the simplest extension of the Hubbard 
model. In this paper we show that the charge-spin fac- 
torized wave function is an exact eigenfunction of the 
Emery model in the strong coupling limit, however the 
spinless fermions ('charges') represent complex objects 
which are the one dimensional analogs of the Zhang- Rice 
singlets. Furthermore, we demonstrate that the Emery 
model can naturally explain the reduction ofpSpectral 
weight for small momenta seen experimentallyl3~li3 and 
also describes some higher energy features of the photoe- 



mission spectra. 

Since the formation of Zhang-Rice singlets is an es- 
sentially strong coupling phenomenon, the approach we 
present in the paper is the pipst suitable method to apply. 
Weak coupling approachestS are inappropriate to capture 
the additional physics due to the presence of additional 
bands. 

The outline of the paper is as follows: In Sec. II we in- 
troduce the Emery model and the canonical transforma- 
tion leading to the strong coupling effective model. The 
spectral functions within this effective model are calcu- 
lated in Sec. Ill, while in Sec. IV the density-density 
correlations relevant to the EELS experiments are dis- 
cussed. 



II. THE EMERY MODEL 

A. Definition of the model 

The Emery model is given by the Hamiltonian Ti = 
T + U + V. For the kinetic part we take the usual tight- 
binding form. 
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where d] ^ and pj^g ^ denote the hole creation opera- 
tors on copper d and oxygen p orbitals at sites i and 
i + S, respectively. The Cu-Cu distance is taken to be 
unity, i are integers and S = ±1/2. The phase fac- 
tors in the hybridization coming from the symmetry of 
the Cu and O orbitals are absorbed in the definition 
of the d and p operators as dj = (— l)-'dphys,j5 and 
Pj+i/2 = (-l)-'Pphys,j+i/2, where the subscript 'phys' de- 
notes the operators respecting the phase factors in the 
hybridization. The inclusion of the phase factors causes 
a shift of TT in the momentum of the O hole and will be ex- 
plicitly mentioned when necessary. In the potential part 
we include the on-site energy difference A = Sp — Sd and 
the on-site Coulomb repulsion U of the Cu 3d-orbitals: 
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C^E<T' 
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where nf,^ = aI,o-a,,o- (« ^ d,,p), and nf = J2a 
thermore the nearest neighbor Cu-0 repulsion 



(2) 



fur- 



(3) 



which may lead to exciton formationtJ. We choose 
L'^ > A in order to haiaj a charge transfer insulator with 
one hole per unit celO. 

Note that if one begins with a model which also in- 
cludes oxygen orbitals on the side of the chain, as a pre- 
liminary step one may build bonding, and anti- and non- 
bonding combinations of these. Therefore, if one under- 
stands our single oxygen orbital per cell to correspond to 



the bonding combination, our results may also be seen to 
represent a good approximation to paiptipfjthe spectrum 
of the more complete four band modelEjO. 



B. Effective model in the strong coupling limit 

As mentioned earlier, direct numerical methods, such 
as exact diagonalization, work only for rather small sys- 
tem sizes. However, in the strong coupling limit (C/, A 3> 
t, V) it is possible to do controlled calculations both an- 
alytically and numerically. As a first step, we derive an 
effective strong coupling Hamiltonian. In the extreme 
case when t = V = 0, the Hamiltonian is block-diagonal 
in the subspace of states with a given eigenvalue of U. 
The hybridization in T may be iijeated perturbatively us- 
ing a canonical transformationcj, leading to an effective 
Hamiltonian TCeS acting within one subspace. A detailed 
and systematic explanation for the case of the Hubbard 
model is given in Rcf. pSl We denote by tilde (e.g. p) the 
operators acting in the subspace of Ticff , and the physical 
operators O are then obtained from 



where O = 0{p -^ p,d ^ d, . . .) and 



(4) 



s^^{r^-Ti 



^—r{'fu^A-'r^-A) + o{t') 



U~A' 



(5) 



is the generator of the canonical transformation with 
TA^-tY^pUsAA^-nU), 



i.S.i 



Trr 



^^E4. 



U-A = -r 2^ ^aPi+S^a^t,. 
i.S.a 



(6) 



The subscript nU+mA denotes that the state acted upon 
is promoted to a subspace at this energy difference. In 
other words 



[U, OnU+mA] = (nU + mA)OnU+r. 



•.A, 



(7) 



and every operator can be decomposed as O = 
X^T! m ^nU+rnA with 11 integers and m integers or half- 
odd integers. Then Tics is given by: 



n,« = u + v + ^ 

+0{t^). 



-'A' -'a 



1 



U-A 



-'[/-A' ^U- 



(8) 



Separating the different processes, Hcfr — W + Ho + 
Til + 7^2, and introducing the effective hopping ampli- 
tudes ts — t^ /{U — A) and It = t^ / A, furthermore 
U' = U + Atr + 4is and A' = A + 4iT, we get (see 
also Refs. M^: 



A' 



i i i,S 

1.(7,(5 

Hi = (ts + It) Y1 {p\+S.JiA,aP^+5' ,a " Pl+S^a^sP^ 

1,8,5' , a 

'H2^tsY^ nf+2S,ad\+2SMK'yKs- 



■i + 6',0" 



^rE^'Wp»-«,<T' 



(9) 



Tio and 7^2 describe the motion of the empty and doubly 
occupied site, respectively, while TLi is responsible for the 
dynamics of the hole on oxygen. By using the Heisenberg 
model ground state |GS) for the insulating|-efise, the main 
effect of the (fourth order) AF interactioncj between Cu 
spins is accounted for. 



III. SPECTRAL FUNCTIONS 

Now let us turn to the spectral functions. The photoe- 
mission spectrum is proportional to the single particle 
spectral function, defined by 

B{k,u) = 5] |(/bLlGS)P<5(c. + Ef - Eos) 

+ Y.Mdl.\^S)\H{Lo + Ef~ Eos), (10) 

assuming that the cross sections of the Cu and O electron 
removal are equal. The sum is over final states |/) with a 
hole added, and a similar definition holds for the inverse 
photoemission spectra, where a hole is removed. First, 
we will present the analytical and numerical calculation 
of the spectral function for the effective model, and then 
we will compare our results to the spectral function of 
t-J model and the photoemission spectra of SrCu02 and 
Sr2Cu03. 



(11) 



PU.;A/2=PW + 0(t'), 

4<t;A/2 = ~X E[(1 ~ f'U)pl+S,a - dlaPl+5,A,s] 
6 

rAJ2^^tfPl+S,a + dlaPl+S,A,a] 



U-A 
+0{t^). 



(12) 



The final states in peaks 'b' and 'c' are obtained by ap- 
plying c[. ^.^ /2 to the ground state, and the sum rule of 
the 'b' + 'c' peaks in Fig. || is 



(GSbfc^,.A/2PL;A/2|GS> 
(GSMfc,o-:A/2'^fc,<T;A/2l'-'S) 



l + 0(i2). 



(13) 



Similarly, the weight in peak 'a' and 'd' is l/2+0{t'^). For 
simplicity, we do not include the corrections in Eq. (|l2| ) 
when we calculate the spectral functions. 
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FIG. 2. Schematic distribution of the weights in the strong 
coupling limit for a charge transfer insulator. 



A. Calculation of spectral functions 

Since the nearest neighbor Coulomb repulsion leads 
only to a uniform shift of the final state energies in 
the spectral functions, it will be neglected in this sec- 
tion. The strong coupling behavior of the photoemis- 
sion spectra is schematically shown in Fig. g: the hole 
can go either to the Cu site (peak 'a') or to the O 
site (peaks 'b' and 'c'). The creation operator c]. ^ — 

dka-iPka entering the calculation of the spectral func- 
tion Eq. (nfl) can be decomposed in leading order as 



-^k.ij 



-/c,o-;-A/2 



-'/c,cr;A/2 



t 
-fe,(T;;7-A/2 



which represent 
a hybridized mixture of Cu and O atomic states. For 
example, including 0{t) corrections 



The dynamics in the peak 'a' is governed by 'H2'. the 
extra hole on Cu, created by d\ ^ in Eq. (|l^), hops to 
neighboring Cu sites with amplitude ig, leaving the spin 
sequence unchanged: 



^2 



^i-i 



2cr,- 



= — isl • • •2crj-icrj 



M 



^j-i^j 



C7,2- 



(14) 



Here '2' in the wave function denotes the position of the 
extra hole (site occupied by two holes) , and aj the spins 
of the singly occupied Cu sites with j — 1, . . . ,L — 1, a,s 
there are L — 1 spins remaining. The situation is identi- 
cal to the case of U/t -^ -l-oo Hubbard model, where we 
know that the wave function of a state, wth momentum 
k factorizes into charge and spin parts 



L-1 



\fQik)) = 4f E e'^^'-^^^IV',) ® \XL-i{Q,nQ)), (15) 



VI 



i=o 



where {ipj) describes free spinless fermions on L lat- 
tice points with an empty site j, which in our case 
is the site with the extra hole. \xL-i{Q,nQ)) is the 
squeezed spin wave function of i — 1 spins with momen- 
tum Q = 2ttJ/{L — 1), J integer, and other quantum 
numbers nq . The energy of the state is 



eQ{k) = [/' - A72 - 2ts cos(fc - Q). 



(16) 



Now that we have both the energy and the wave function 
of the final state, we are ready to calculate the spectral 
function as presented by Sorella and Parola for the large- 
U Hubbard modelQ. As a first step, we write the ground 
state also in a product form 



|GS) = |Vgs)®Ixgs), 



(17) 



where the |xgs) is the Heisenberg ground-state wave 
function and \ipGs) is the fully filled Fermi see of spinless 
fermions (charges). It is convenient to choose systems 
with 2,6,10 etc. sites, where the momentum of the |GS) 
is zero. In the matrix element of Eq. l fiX\ ) it suffices to 
keep the momentum dependence in the final wave func- 
tion only. 



\{fQik)\dljGS)\' ^ L\{fQ{k)\d]^,JGS)\' 

= |(x(Q,"-q)I^o.<t|xgs)| 



(18) 



where we have substituted the factorized form Eq. ( p^ ) 
and used the fact that the overlap in the charge part is 
1. Only the spin part is nontrivial: the operator Zj_cr 
removes a spin a at site j, reducing the spin sequence to 
length L — 1. Introducing 



DiQ)=^J2\^xiQ,nQ)\Zo^a 



ixGs; 



"Q 



for the spectral function we get: 

i?(fc,c.) = ^i?(Q)5(c. + eQ(fc)). 



(19) 



(20) 



D{Q) is essentially the "occupation number" of the 
spinous, and has a singularity at the spinon fermi mo- 
menta Q = ±7r/2. It can be approximated as (L — 
1)D{Q) « -0.5 + 111^ ^7r2 - 4Q2 for -7r/2 <Q <tt/2 
and zero otherwisa3'QE3. We therefore find that the spec- 
tral function in the upper Hubbard band (peak 'a') is 
identical to that of the large-[/ Hubbard or small J t-J 
model. Note also that the inverse photoemission spec- 
trum will have a similar form with bandwidth Atx (peak 
'd' in Fig. I). 

Let us now proceed to peaks 'b' and 'c' in Fig. 0, which 
can be associated with the hole on oxygens. The hole 
added to an O site with pj, ^ can form a singlet or triplet 



with a neighboring Cu spin. We will denote by \S) {\S)) 
states where the O hole forms a singlet with the Cu spin 
on its right (left), as seen in Fig. || For example, 
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FIG. 3. The state | i 5 Ti>- The Cu and O spins in the 
shaded region form a singlet, which we denoted by S. 

As we will see below, a suitable combiipiion of these 
states will give us the Zhang-Rice singletal3, in terms of 
which the lowest energy excitations may be described by 
a one-band model. In the present calculation we are also 
including the triplets and high-energy singlets in order to 
describe higher energy excitations. Note that this basis 
is not orthogonal: 



{■■■Sa---\---Sa---)^l, 
{■■■Sa---\---(jS---) = 1/2. 



(21) 



In general, the resulting spectrum is complicated and the 
singlets and triplets mix with one oth erj ,, e xcept for the 
particular case of ts finite and iy — CcjS Due to the 
very special form of Tii for iy = 



Hi = -2U 



i,5,S' ,a^a' 






Tga' qri \ ^^f 



i j Pi+S,aPi+S',a'' 



(22) 



where we can identify the projector onto spin singlets 

(here Sf = X^acr' 'A a''''^<^"^i a'/'^ ^"^^ ■'" i^ ^^^ vector of 
Pauli matrices), the matrix elements of Tieff leading to 
propagation of the singlets are: 

ni\---Sa---)^ts\---a(S-^)---) 

-2ts\---{S~*S)<7---), 

+2ts\---<j{^-^)---). 

The combination | S) — \ S) moves through the lattice like 
the site '2' in Eq. (|l|) 

Hil • • •a,_i(^- S)a, • • •) = ts(l • • •c7,-iaj(^-^) • • •) 

+ I-- -(5 -S)a,_ia, •••)), 
(23) 



leaving the spin sequence unchanged. In this case ji/'j) in 
Eq. ( [l5|) will denote this particular combination at site 
j. The energy of the state is 



eg(fc) = A/2 - 4^5 + 2ts cos(fc - Q). 



(24) 



Next, we need to calculate the matrix elements. Using 
the identity 

|(/|pl.,JGS)p=i|(GS|pi/2,.|/(fc))P, (25) 



where ^1/2.0- removes the hole at site i = 1/2, the k de- 
These singlets leads to the formation of the 'c' peak in pendence is now in the final state only. So, for the matrix 



Fig. |. 



element we get: 



j 
= (1 + cosk)\(xGs\zl^\xiQ,nQ))\'' , 



(26) 



where we have used that (xcsl^J tlxlQ)) = 

—e^'^{XGs\Zl-t\xiQ))- Dividing the matrix element by 
the norm of the final state 

A — 2e^ (k) 
(/Q(fc)l/Q(fc))=2-cos(fc-g) = ^^, (27) 

and summing over final states with definite Q, we can 
write the spectral function as 



Bsik,Lu) = 



4is(l + cosfc) 

A + 2uj 



Y,DiQ)S{Lu + el^{k)). (28) 



Clearly, even introducing form factors in the one-band 
model [which is identical in form to Eq. (po|)], the w- 
dependent prefactor of the spectral distribution above 
cannot be obtained. The local (fc-averaged) spectral 
function for the singlets is 

^ ^ 1 4tg + (21n2-l)(-2c^-A+8^g) 

^ TT {A + 2uj)^{A~Ats + 2uj){~2uj-A + 12ts) ' 

(29) 

with weight 0.32. The rest of the weight (0.68) is at 
higher energies uj — —A/2, where wc find nondispers- 
ing solutions, made of a particular combination of sin- 
glets \Saa') - 2\(7Sa') - 2|cr5cr') + {aa'S), as weh as the 
triplets, contributing with a delta peak to the spectral 
function to form the peak 'b' in Fig. 0. 

The only requirement for the procedure outlined above 
to work is that during the motion of the hole the spin se- 
quence is unchanged. This immediately requires tx = 
in Hi- 

The influence of finite tr is shown in Fig. ^. The lower 
'singlet' band increases its width, while the overall shape 
of Bs{k,uj) does not change significantly. On the other 
hand, peak 'b' now extends from A/2 to A/2 + 4tT and a 
sharp dispersion dominates the spectrum. Only a slight 
weight transfer from the 'singlet' to the 'triplet' band can 
be observed, e.g. a,t k — {q — tt) the weight in the 'sin- 
glet' band is reduced from 0.65 to 0.43 as we increase t^ 



from to ts, while the total weight, given by Eq. (|l^) is 
unchanged in leading order. 



B. Comparison with the t-J model and 
photoemission experiments 

Comparing with the t-J model for small J (Fig. |p,), 
we can see that although (also for t^ = ift) the 'sin- 
glet' feature is similar to the t-J model resullfi, there are 
detailed differences in the distribution of weight, simi- 
lar to those in Eq. (p8|), as well as in the dispersion of 
the upper edge of the 'singlet' continuum. We therefore 
see that even in parameter regimes where the one-band 
t-J description accurately predicts low-energy excitation 
energies, the two-band model may have significantly dif- 
ferent properties as far as other physical observables is 
concerned, exemplified here by the momentum and fre- 
quency dependence of the spectral weights. The effect 
of finite J is to give dispersion to the now dispersionless 
lower 'spinon' edge in both Emery and t-J model. 

Finally, let us compare our results with the experi- 
ments. For both SrCu02 and Sr2Cu03 the low energy re- 
gion shows features found in the t-J model, i.e. the holon 
and spinon bands dispersing with t ss 0.5 — 0.6eV and 
J « 0.15 — 0.2eVp. respectiicely, which is consistent with, 
the susceptibility^, opticaO and electron-energy lossllll 
experiments. However, an additional interesting feature 
is the weight reduction as the zone center {q = 0) is ap- 
proached. In Refs. [7|,p^ this is attributed to the different 
cross sections of Cu and O orbitals, while in our theory 
it arises quite naturally from the internal structure of 
the low energy singlets. Concerning the higher energy 
features, the 'triplet' feature is in reasonable agreement 
with the dispersing peak at 2 eV in Fig. 6 of Ref. pi if 
one disregards the flat non-bonding oxygen bands not in- 
cluded in our model. These general trends don't depend 
strongly on tx/ts and the inclusion of t/A and t/{U — A) 
correction in operators leads only to a small weight trans- 
fer to lower energies. 
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IV. DENSITY-DENSITY CORRELATIONS 



This leads to the sum rule 



The density-density correlation function describes the 
dynamical dielectric response of the material and is acces- 
sible by measuring e.g. the optical conductivity, electron- 
energy loss spectra (EELS), aiad inelastic X^jay scatter- 
ing. Both optical conductivityEJ and EELSEJ have been 
measured on Sr2Cu03. The EELS spectra can be reason- 
ably well interpreted within a one band Hubbard model 
extended with nearest neighbor repulsion, so it is inter- 
esting to see what changes if we consider a two band 
model, like the Emery model. 

Since the system is an insulator, we get finite density 
response only above the charge transfer gap at w « A. 
In lowest order the density-density correlation function 
is given by: 

J^{k,uj) = ^ |(/|nfe;A|GS)p5(c^ -Ef + Egs), (30) 



^2 ^ 

7V(fc, Lu)duj = 8-^ sin^ -. 



(32) 



Now let us determine the dynamical density response 
in the special case tx — 0. The operator nk;A moves a 
hole from Cu to O, and results in a two body problem, 
which can be solved using standard techniques. For the 
'singlet' part, the final state wave function can be repre- 
sented as 

L-l 



j=i 



with the norm 



(33) 



where rik-A = ^t-A + "^fe-A '^^^ ^^ calculated from rii-A 
(l/A)[7A,n.j] and reads 
t 1 



nk-A 






,:St 



A/L 



e-^-l)(l-n,-,)p]+,^,.^.,. 



J,S,cr 



(31) 



{fs\fs) = J2(^ 



L-l L-2 

y'j)-YyjXj+i- 



(34) 



These states are L-fold degenerate, since 'e', which rep- 
resents the Cu with no hole, does not hop in this limit. 
The Schrodinger equation gives 



Exj = Eyj == ts{xj+i - 2x.j - 2y.j + y^^i) (35) 

for j = 2,- ■ ■ , L — 2 with boundary conditions 

Eyi = ts{x2 - 2xx - 2yi), 
E{yi - xi) = Vxi, 

EyL^l = ts{xL-2 - 2XL-1 - 22/L-i), 

E{yL-i - XL-i) = VxL-i- 

The energy E is measured from A. Due to symmetry 
Xj and yj are real and there are even and odd parity 
solutions with x. 



0. We immediately 
: yj for £' 7^ and 



-J - TVL-j- 
Let us first consider the case V 
see that (xj — yj)E = 0, i.e. Xj 
J = 1, • • • , (i — 1). The solution is Xj = yj = sin jk, with 
K = It^/L, and / = 1, 2, • • • , L — 1. These states have 
even (odd) parity for / even (odd), energy 



E,. 



~Ats + 2ts cos K, 



(36) 



and norm {fs\fs) = L{2 — cosk)/2. For L ^ oo they will 
form a continuum from E = —6ts to — 2<5. Additionally 
there are L—1 degenerate states with E ^ 0. The matrix 
elements in Eq. ( |30| ) for the singlet contribution read: 



oy-2 
|(/sK;A|GS)P = ^F5 






1 ± COS ■ 



(37) 



where the + (— ) sign is for the even (odd) state, and 
Fs = (GS|i - So • Si|GS) -^ ln2 for L ^ oo is the 
probability of finding two neighboring spins forming a 
singlet in the spin sequence. The only nontrivial quan- 
tity is Xi/{fs\fs), which can be conveniently expressed 
using the energy of the state as 



(/l/> 



1 2 sin K 



L2 



1 {2ts + E) {6ts + E) 
L tsE 



(38) 



In the thermodynamic limit we replace the sum over 
states with an integral over energy: 



E 



dE^L''^ 



'TTdEj 



(39) 



where {l/'K){dK/dEK) is the density of states 
1 Sk 1 1 



Tra^K t: ^{^2ts-E){ms+E)'' 



(40) 



the factor [1 ± cos(A:/2)] in Eq. ( pTJ ) averages to 1, and for 
the contribution of the Zhang- Rice singlets to the density 
response we get: 



,^,, , 81n2 i^ ^(u;-A + 6is)(A-2ts~w) . 2^ 
^^k,u;) = --- ^^^— ^ sm -. 

(41) 

The density response in this limit has a trivial momen- 
tum dependence, due to the nondispersing nature of the 
'e' site. It gives 2(2 - \/3)Fs « 37% of the total weight, 
the rest of the weight is in a single peak at lj = A, which 
comes from the nondispersing singlets and triplets. The 
A/'(fc, uj) we just calculated is shown on the upper left plot 
in Fig. I 

Turning on the Cu-0 repulsion, which acts as an ef- 
fective attraction between the empty Cu site and the O 
hole, two twofold degenerate (for L — > oo) excitons with 
energies 



"1 = 



^tsV 



V^±Vy/l2tl 



v^ 



ts-2V 



(42) 



appear, together with a twofold degenerate exciton in- 
volving the triplets at ^t — —V. The J7g solution exists 
only for V > 2ts where it splits off from the lower edge of 
the continuum. Not going into the details, the expression 
for x\/{fs\fs) in Eq. (|37|) now reads 



xl 



-E {2ts + E) (Gts + E) 



iflf) 24t|y + 2E{Ats - V)V + L [{2V - ts)E^ + 2(3ts + V)VE + 3VHs 



(43) 



which is valid both for the f2g excitons and the contin- 
uum. The Nik, to) is complicated and we do not give 
the analytical form, which is straightforward to derive 
from Eqs. (|^), (Q, and (H), but we refer to Fig. | 
for a discussion of features. For small V the energy of 
the exciton is f2g w (3 — 2^/2,)V with relative weight 
« [2V3 - 3 - 2(7 - 4v%)y/is] Fs, i.e. increasing V the 
weight is transferred to the continuum. For large repul- 
sion {ts ^ V <^ U) the weight is concentrated in the 
—V— 2ts, while the continuum and 



excitons at ilj 



the exciton at fi^ ? 
the order of t|/V^^- 
citon fir has weight 1 
size of V. 



y —'its/2 has a negligible weight of 
On the other hand, the triplet ex- 



Fs ~ 31%, independent of the 



To study the effect of finite It, we used numerical 
Lanczos diagonalization of small systems (18 site) to ex- 
tract the density-density correlations. As can be seen in 
Fig. ^, the effect of finite ir is substantial: (i) Because 
'e' acquires dispersion, the shape of the spectrum resem- 
bles more closely tiiat of the one-band model, in that it 
narrows as fc — > tm. Turning on V , excitonic features 
are formed at the zone boundary, (ii) The contributions 
coming from 'singlets' and 'triplets' cannot be separated 
energetically, especially for small momenta, (iii) While 
for ^T = the excitons are sharp peaks, for ty 7^ these 
sharp excitonic peaks broaden and form an incoherent 
spectrum. A direct comparison with the experiments to 
decide whether a two band model is more appropriate is 



difficult, given that the resolution of the EELS spectra 
is not sufficient to see the detailed differences between 
the two band and the one band models. Furthermore, all 



of our results here are valid in the strong-coupling limit, 
so one expects quantitative changes to the spectra for 
experimentally relevant coupling strengths. 
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FIG. 5. Nik, uj) of 18 site effective model for V = G, 2ts and 8fs from left to the right, and tr = G (upper) and tr ~ ts 
(lower plots) obtained by exact diagonalization. The thick line in the upper plots shows the analytical result. The S functions 
are plotted as Lorentzians of width O.lts and the plots for each k are normalized to have total weight 1. In the upper plots 
we show k — TT only because of the trivial k dependence. The remnants of excitons are indicated within the parenthesis in the 
lower right plot. These spectra are independent of the phase factors of Cu and O orbitals. 



Recently, the density response was calculated using 
projection techniques in Ref. p2l with a result which dis- 
agrees with ours. Since the calculation presented in that 
paper is rather involved, it is difficult to trace whether 
the difference is due to the method applied or to the pa- 
rameter regime investigated. 



momentum dependent corrections. This provides a very 
simple and natural explanation for the momentum and 
frequency dependence of the spectral weights observed 
experimentally. 



V. CONCLUSION 
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